NIC  FILE  COPY 


LEVEt^ 


EFFECT  OF  RECORD  LENGTH  ON  THE 


CORRELATION  OF  COMPLEX  EXPONENTIALS 


Tapan  K.  Sarkar 
2 

Joshua  Nebat 

2 

Donald  D.  Weiner 
VI jay  K.  Jain3 


(  # 


This  work  has  been  supported  by  the  Office  of  Naval  Research  under 
contract  N00014-79-C-0598 . 


1.  Dept,  of  Electrical  Eng. 

Rochester  Inst-,  of  Technology 
Rochester,  New  fork,  14623 


3.  Dept,  of  Electrical  Eng. 
University  of  South  Florida 
Tampa,  Florida  33620 


2.  Dept,  of  Electrical  Engineering 
Syracuse  University 
Syracuse,  New  York  13210 


Approved  io:  public  velvoim; 
Diatribe  tic  >a  U  nH/riled 


*1  11 


in 


SECufclTY  CLASSIFICATION  OF  THIS  PAGE  fHTim  Dote  Enlarad) 


REPORT  DOCUMENTATION  PAGE 


.  REPORT  NUMBER 

TR-81-4 


4.  TITLE  (end  Submit) 

Effect  of  Record  Length  on  the  Correlation 
of  Complex  Exponentials 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


J.  RECIPIENT'S  CATALOG  NUMBER 


3  TYPE  OF  REPORT  B  PERIOD  COVERED 

Technical  Report  No.  4 


performing  org.  report  number 


t.  author^; 


J.Nebat,  T.K.Sarkar,  D.D. Weiner  and  V.K.Jain 


».  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Dept,  of  Electrical  Engineering 
Rochester  Institute  of  Technology 
Rochester,  New  York  14623 


II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 


Department  of  the  Navy 
Office  of  Naval  Research 
Arlington,  Virginia  22217 


4.  MONITORING  AGENCY  NAME  A  AODRESSflf  dllltrent  tram  Controlling  Olttco) 


•  •  CONTRACT  OR  GRANT  NUMBCRfaJ 


N00014-79-C-0598 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 

61153N  RR021-01 

RR021-01-01 

NR371-014 


*2.  REPORT  OATB 

November  1981 


is.  NUMBER  OF  PAGES 

21 


IS.  SECURITY  CLASS,  (ol  tbit  report) 


I 


UNCLASSIFIED 


a.  oeclassification/downgraoing 

SCHEDULE 


17.  DISTRIBUTION  STATEMENT  (ol  a It  tbttrtcl  tnttrtd  In  Block  30,  II  dllltrent  bom  Report) 


IS.  KEY  WORDS  (Continue  on  rtvtrtt  tldt  II  neceeeory  tnd  Identify  by  block  ntmbor) 

1)  Correlation 

2)  Decaying  Exponentials 

3)  Orthoganal ity  of  Exponentials 

4)  Resolution  of  Exponentials 

5)  Time  Domain  Identification 


20.  ABSTRACT  fConffnua  on  ravaraa  old*  if  nac»a««ry  and  identity  by  block  number) 

Resolution  of  two  sinusoidal  signals  of  frequencies  f^  and  f^  requires  a  time 
record  of  length  at  least  1/ | f ^-f ^ | -  Resolution  of  two  decaying  exponentials 
of  complex  frequencies  Sj=a^+juj^  and  with  correlation 

coefficient  is  dependent  only  on  the  ratios  c^/o^,  and  for  a  time 

record  of  length  greater  than  2f\a^\.  This  is  also  the  condition  for  near- 


DD  ,  jan^ts  1473 


COITION  OF  I  NOV  SS  IS  OBSOLETE 
S/N  0102-LF-014-4601 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  TNIS  P  AGE  fWRan  Dot o  knitted) 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  FAOt  (*htn  D«(  Bnfnd) 


20. 

orthagonal ization  of  a  set  of  complex  exponentials,  with  small  error. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  FAOin»Mn  Omm 


1.  ^Introduction 

With  the  advent  of  the  Singularity  Expansion  Method  (SEM)  there 
has  been  a  great  interest  on  the  identification  of  a  linear  time- 
invariant  system  by  a  sum  of  complex  exponentials.  In  this  paper  the 
suitability  of  the  exponential  functions  for  modelling  a  finite  time 
domain  impulse  response  is  examined.  More  specifically,  we  address 
the  question  of  how  long  a  data  set  does  one  need  so  that  the  record 
length  behaves  as  if  it  were  infinite.  In  other  words,  what  is  the 
minimum  length  of  record  required  to  resolve  the  various  components  of 
decaying  exponentials.  The  answer  to  this  question  may  yield  data  for 
a  meaningful  analysis. ^ 

If  two  signals  <(>(t)  and  ¥(t)  are  to  be  distinguishable  the  waveforms 
must  have  the  property  of  being  as  different  from  its  shifted  self  as 
possible.  In  mathematical  terms,  the  mean  squared  departure  of  YCt) 
from  <J>(t+r) 

CO 

/|nt)-$(t+T)j2dt  (i) 

— oo 

must  be  as  large  as  possible  over  the  range  of  T.  By  expanding  the 
above  integral  and  noting  the  independence  of  the  squared  terms  of  t, 

we  see  that  minimization  of  (1)  implies 

00 

fV  (t)4>*(t+r)dt  (2) 

_oo 

shall  be  as  small  as  possible.  Here  *  denotes  complex  conjugate. 

The  above  integral  in  (2)  is  defined  as  the  correlation  between 
the  two  functions  'i'(t)  and  4>(t).  We  now  introduce  a  normalized  version 
of  (2)  which  we  define  as  the  correlation  coefficient  p(x)  between  the 


two  signals  4>(t)  and  'F(t).  The  correlation  coefficient  p  (x)  is  defined  as 


2 


CO 


Jv(tw*(t+T)dt 

Pjx)  =  - - If - - - - - - -  .  (3) 

{  /H'(t)4'*(t)dt}'2  {  /4>(t)4>*(t)dt}^ 

.00  .00 

So  for  perfect  correlation,  i.e.  when  the  waveforms  VCt)  and  <j>(t) 
are  identical  then  p  (t)“1.  Under  this  circumstance  it  is  impossible 
to  resolve  H'(t)  from  4>(t).  However,  when  p  (t)*=0,  we  have  perfect 
resolvability,  i.e.  the  two  signals  ^Ct)  and  $(t)  are  as  different  as 
possible. 

Observe  that  the  limits  in  the  integrals  of  (3)  are  from  to 
+<».  In  general  when  we  are  performing  an  experiment  it  is  not  possible 
to  have  infinitely  long  data  records.  If  we  have  finite  length  data 

records  then  the  correlation  coefficient  is  defined  as 

T+A  * 

/ntH  (t+T)dt 

PA(T)  =  T+A  ~  .  T+A  “  T  ' 

{  /vct)^  (t)dtP  {  /  4>(t)<f>  (t)dtP 

T  T 

It  is  clear  that  because  of  finite  record  lengths, 

|pA<t>I  f  IPoo^I’  (5) 

In  this  paper  we  investigate  the  value  of  A  for  which 

!  PA(t>  I  =  I  Poo^T>  I  •  <6> 

This  value  of  A  will  then  dictate  the  length  of  record  necessary  to  reduce 
the  correlation  coefficient  between  the  functions  'F(t)  and  <f>(t)  and 
thereby  increase  the  chances  of  resolvability. 

2.  Correlation  Considerations  Involving  Complex  Exponentials 

In  general,  identification  of  the  complex  exponential  components 
of  a  signal  involves  solution  of  a  set  of  simultaneous  equations  (e.g. 
Prony's  method).  As  the  correlation  between  the  components  is  increased, 
the  equations  become  more  ill-conditioned.  Consequently,  the  correlation 


3 


coefficient  between  the  two  complex  exponentials  is  a  measure  of  the 
difficulty  by  which  the  two  components  may  be  resolved.  By  way  of 
illustration,  two  signals  with  unit  correlation  result  in  a  singular 
set  of  equations  whereas  two  signals  with  zero  correlation  yield  an 
uncoupled  set. 

Consider  two  simple  complex  exponentials  given  by 


{eSit};  i  =  1,  2, 


0  <  t  <  00 


where  Sj*<Jj+jui  ,  with  0  <0  and  j=  /^T.  The  correlation  coefficient 
between  exp  [s^t]  and  exp  [s2t]  over  the  time  interval  [T,  T  +  A] 


is  defined  as 


T+A  s.t  s  *t 

t  1  2. 

j  e  e  dt 
T 


Pa(t) 


A  s.t  s  *t  ,  T+A  st  s  *t  , 

{  /  e  1  e  1  dtP  {  /  e  2  e  2  dtP 


Performing  the  integration  and  rearranging  terms,  we  get 


pa<t) 


e(CT2^2)T.eJS~a)2)T-{e[(0l'KT2)+j(aira)2)lA_1} 


[Oj+t^+j 


In  general  the  correlation  coefficient  is  a  complex  quantity.  The 
initial  time  instant  T  simply  adds  a  constant  to  the  phase.  Since  we 
are  primarily  interested  in  the  magnitude  there  is  no  loss  in  generality 
by  assuming  T=0.  Also  the  correlation  coefficient  varies  as  exp{  [o-ju^lx} 
with  t,  which  does  not  enter  into  our  discussions.  Hence  we  define  a 
new  coefficient  Pj2(A)  as 

PA<x)  -  p12(A)exp{  [o2-ju)2]r}  (10) 


We  want  to  study  the  properties  of  Pj2(A)  as  the  exponential  function 
in  (10)  does  not  provide  any  additional  insight. 


1 


Of  particular  interest  is  the  sinusoidal  case  for  which  0^=0  for 
i=l,2.  In  this  case 


sinIW]/(^2,  A 
P12(A)lo1  2=0  “  (oj1-u2)A  j  2 


The  phase  and  the  magnitude  of  the  correlation  coefficient  is  plotted 
in  figure  1  for  to  Observe  that  the  correlation  coefficient 

is  zero  for 


(o^-oo^A  2w(f1-f2)A 


=  nTT ;  n=l  ,2, . . . 


Hence,  for  an  observation  interval  of  length  A,  the  two  components 
are  uncorrelated  for 

frf2=  ?  5  n  =  1»2*-”  •  (13) 

With  n  =  1,  we  obtain  the  conventional  condition  for  resolution  of  two 
frequency  components  which  is  given  by 

frf2-  i  •  <“> 

for  other  values  of  fj-f2,  t*ie  magnitude  of  the  correlation  coefficient 
is  bounded  by 

|p  CA)|  0  1  i  (15) 

12  Va2  TT(f1-f2)A 

This  bound  is  plotted  in  Fig.  1.  When  a)j=aj  ,  the  correlation  coefficient 
is  unity.  This  is  to  be  expected  since  the  two  components  are  then 
identical. 

The  situation  is  more  complicated  when  <  0;  1=1,2.  This  is 
evident  by  consideration  of  the  spectra  involved.  For  the  infinite 
interval,  0^=0  yields  a  line  spectrum  whereas  o^<0  results  in  a 
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continuous  spectrum.  Define 


a  ~  a  ;  ® 
°1 


“l""2 

—  5  Y  "  °1A 


where,  for  convenience,  it  is  assumed  ° 2  —  a i  ~  Substitution  of 
(16)  into  (8),  with  T=0,  the  correlation  coefficient  is  expressed  as 

r—  (l+a)Y  jSY  . 

Pl-(A)  ~  fmnc  - - - - —  07) 

12  1+ot+jp  ^  2aY^ 

(l-e^Y)(l-e^aY) 

We  define 

| p12 (A) |  =  |p12(«)| *W(a,8,Y)  (18) 

where  p^C00)  the  correlation  coefficient  due  to  the  infinite 
observation  interval  and  is  defined  as 


Ip12(°°) 


I  2  2 

(1+a)  +6 


The  second  factor  in  (18)  can  then  be  interpreted  as  a  "windowing” 
factor  due  to  finite  observation  interval.  The  window  factor  is  defined  as 


W(a,8,y) 


i_(e2V0‘W<i'K‘)Y 


For  a  fixed  value  of  $,  the  maximum  value  of  occurs  for 


ctn  -'Vi+e 


The  peak  value  of  |  p ^ 2  C°°)  I  *s  then  given  by 


lp12(0°)ia=aM 


1  +  ^1+8' 


It  is  interesting  to  note  that  (17)  is  unchanged  by  interchanging  the 
subscripts  1  and  2  in  the  definitions  of  (10).  For  this  reason,  it 
is  necessary  to  consider  only  values  of  a  greater  than  or  equal  to 
unity.  The  value  of  a  -  1  is,  therefore,  of  special  interest  because 
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it  serves  as  the  "origin"  in  our  plots.  Interestingly  enough,  for 
both  large  and  small  values  of  |B|,  the  peak  value  in  (22)  is 
approximately  related  to  the  value  of  the  correlation  coefficient  for 
a  =  1.  Specifically, 

4 


lp12(00)la=l 

for  1 6  j  »  1, 


|p12(">la.l  “W” 


On  the  other  hand,  for  |g|  «  1, 

lp12(”>la.l  “  1  "  lon^L. 


aM 


(23) 


(24) 


(25) 


Also,  for  a  >>  1  and  a  >>  |fj|,  observe  from  (19)  that 

lp12(<X>)l  *  ~k  (26) 

Hence,  for  a  fixed  value  of  |B|,  the  correlation  coefficient  approaches 
zero  in  the  limit  as  a  approaches  infinity. 

Equation  (19)  is  now  investigated  as  a  function  of  8.  For  a  fixed 


value  of  a,  the  maximum  value  of  Ip^C00)!  occurs  for 

6m=°- 


(27) 


It  is  interesting  to  note  that  the  correlation  coefficient  peaks  when 
<i)l=UJ2  kUt  does  not  necessarily  peak  when  0^=0  (i»e*»  ct=l)  •  For 


B„=0,  the  peak  value  of  the  correlation  coefficient  is 
M 


P12(°°)| 


-  1*5 

B=Bm=0  ~  14a 


For  a  »  1,  the  peak  becomes 

2 


|P12<”>I 


Asymptotically,  for  |B|  »  a,  (19)  reduces  to 


|P12(°°)|  = 


l/a. 

TeT 


(28) 


(29) 


(30) 
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Therefore,  for  a  fixed  value  of  a,  the  correlation  coefficient  approaches 
zero  in  the  limit  as  j 3 )  approaches  infinity.  By  comparison  of  (26) 
with  (30),  it  is  seen  that  the  correlation  coefficient  asymptotically 
approaches  zero  at  a  faster  rate  with  respect  to  |B|  as  opposed  to  a. 

A  plot  of  |  p ^ 2 (°°)  l  versus  1 3 1  ,  with  a  as  a  parameter  is  shown  in 
Fig.  2.  Recall  that  the  larger  is  the  correlation  between  signal 
components,  the  more  ill-conditioned  are  the  equations  which  arise  in 
the  identification  problem.  Figure  2  points  out  that  the  problem  of 
resolution  is  eased  under  situations  of  both  large  a*s  and  large  B’s. 
Recognizing  that  a  is  constrained  to  be  greater  than  unity,  small  a 
implies  a  =a  .  Then  a  large  value  of  |B|  is  desirable  so  that  the 
difference  in  and  will  aid  in  discriminating  between  the  two 
components.  On  the  other  hand,  if  j  3 1  is  small,  a  large  value  of  a  is 
desirable.  It  is  interesting  to  note  that,  for  large  a,  the  correlation 
coefficient  is  relatively  insensitive  to  | 3 |  (e.g.,  see  curve  with 
a  *  100) .  This  is  reasonable  since  a  large  value  of  a  implies  that  one 
component  decays  much  more  rapidly  than  the  other.  Hence,  the  correlation 
coefficient  is  more  influenced  by  the  relative  decays  as  opposed  to  the 
relative  oscillations.  In  the  identification  problem  a  and  B  are 
specified  and  the  observation  interval  is  finite.  Since  an  infinite 
observation  interval  was  assumed  in  obtaining  the  curves  in  Fig.  2,  they 
serve  as  a  lower  bound  on  the  correlation  coefficient  for  the  finite 
interval. 

A  second  way  of  viewing  *s  presented  in  Fig.  3  where  the 

magnitude  is  plotted  as  a  function  of  a  with  |B|  as  a  parameter.  The 
conclusions  arrived  at  from  the  previous  figure  are  still  valid. 


8 


However,  Fig.  3  clearly  shows  that  large  values  for  both  a  and  | 6 | 
are  preferable.  Also,  note  the  sequence  of  peaks  at  values  of  ot^ 
predicted  by  (21).  For  large  values  of  |8|>  0^= | B | -  The  maximum 
value  of  this  correlation  coefficient  then  arises  when  [  cr^  [  -  loJj-u^l 
where  it  is  assumed  | |  >  |cJjJ. 

Still  another  way  of  presenting  the  results  is  to  construct  equal 
correlation  magnitude  contours  as  a  function  of  a  and  |B|.  Solution 
of  (19)  for  j  6  j  results  in 


|B|  =  2a  (— . . . j-  -1)  -  (1-Kx2).  (31) 

^  |p12(oo)i 

The  contour  plot  is  shown  in  Fig.  4  where  the  parameter  is  J  p ^ 2 C °°) ]  - 
This  plot  enables  the  user  to  determine  the  various  combinations  of  a 
and  | 6 |  which  result  in  a  given  value  of  the  correlation  coefficient. 

The  plot  also  allows  one  to  determine  the  sensitivity  of  the  signal 

parameters  to  small  changes  in  the  correlation  coefficients.  For  example, 
assume  0^=3  and  ]8|=4.  For  J  (°°)  I  =0. 4 .  the  allowable  value  of 
is  approximately  23,  as  read  from  Fig.  4.  The  corresponding  value 
of  a 2  is  -69.  When  |  p^^00)  l=0*  5,  the  allowable  values  for  a  are  1.3 

and  13.  This  yields  values  for  of  3.9  and  -39,  respectively.  In 

this  case,  the  value  of  is  seen  to  be  highly  sensitive  to  changes 
in  the  correlation  coefficient. 

Our  discussion  thus  far  has  dealt  with  the  infinite  interval.  We 
now  consider  the  effect  of  a  finite  observation.  The  windowing  effect 
is  accounted  for  in  (18)  by  the  second  factor  which  is  W(a,B,'Y)»  and 
can  be  shown  to  approach  unity  for  all  possible  choices  of  a  and  B»  as 
| y |  -*■  °°.  Obviously,  the  effect  of  windowing  is  negligible  when 
W(a,B,Y)-l*  Since  a  1  and  y<0 ,  W(ot,8,y)  is  bounded  by 
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(32) 


<W(a,8,Y)< 


1+e 


2Y 


1-e 


2Y 


Observe  that  the  bounds  are  independent  of  a  and  0.  The  inequality  in 
(32)  can  be  used  to  obtain  an  estimate  for  the  minimum  length  of  the 
observation  interval  in  order  that  the  correlation  between  the  two 
components  be  approximately  the  same  as  for  the  infinite  interval. 

(In  general,  windowing  tends  to  increase  the  correlation).  From  (32) 
the  following  table  is  obtained. 


Table  1.  Lower  and  Upper  Bounds  on  W(a,8,Y)* 


-1 

.930 

1.313 

-2 

.991 

1.037 

-3 

.999 

1.005 

It  is  concluded  from  Table  1  that  the  record  length  can  be  assumed 
to  be  infinite,  as  far  as  W(a,B,Y)  is  concerned,  provided  y  <  -2. 

In  other  words,  when  the  length  of  the  observation  interval  is  such 
that 

A  >  —  where  a .  <o  <0;  i  =  2,3,...,  (33) 

—  i  1 

then  ( p ^ 2 (^>)  |  -  Since  a  finite  interval  tends  to  increase 

the  correlation  coefficient  (see  Figs.  5,  6),  it  is  desirable  that  the 
inequality  in  (33)  be  satisfied. 

The  behavior  of  | p ^ 2 (^) j  as  a  function  of  | Y 1  is  illustrated  in 
Figures  5  and  6.  In  Figure  5,  a  is  constrained  to  be  unity  as  the 
parameter  |0|  is  varied  from  1  to  100.  Note  that  the  curves  have 
settled  down  to  their  asymptotic  behavior  for  ( Y |  >2.  It  is  interesting 
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to  compare  the  damped  case  (i.e.,  0)  to  the  sinusoidal  case 

(i.e.,  0 j  =  0).  For  0^=  -1,  |B|  =  •  The  dashed  line  in  Fig.  5 

corresponds  to  the  sinusoidal  case  where  “  10.  This  serves 

as  a  reference  for  the  damped  case  where  |B|  -  10  and  0^  =  -1.  In 
Figure  6,  | B |  is  constrained  to  equal  the  value  3  as  the  parameter  a 
is  varied  from  1  to  100.  Again  asymptotic  values  have  been  reached 
for  |y|  >  2.  The  dashed  curve  in  Fig.  6  represents  the  sinusoidal  case 
with  | I  =  3.  Since  ( B |  =  3,  each  curve  in  the  figure  may  be 
compared  with  the  sinusoidal  case  provided  is  assumed  equal  to  -1. 
for  a  >  1,  it  is  interesting  to  note  that  there  exist  values  of  |y| 
where  the  correlation  coefficient  is  smaller  than  that  of  the  sinusoidal 
case.  The  only  exception  is  for  a  =  1.  Notice  that  the  largest  asymptote 
occurs  for  a  =  1+B^  as  predicted  by  (21). 

3.  Orthogonality  Property  of  Complex  Exponentials 

In  general,  two  complex  exponentials  as  given  by  (1)  are  not 
orthogonal  over  any  interval  within  [0,°°).  In  this  section,  we  would  like 
to  investigate  under  what  conditions  two  complex  exponentials  (which 
are  orthogonal  in  the  infinite  interval)  remain  almost  orthogonal  in 
the  finite  time  interval. 

The  Gram  Schmidt  orthogonalization  procedure  can  generally  be  used 
to  orthoganalize  a  set  of  functions.  However,  a  simpler  procedure  for 
exponential  functions  was  developed  by  Kautz  [l]  and  applied  in  [2-3]. 

The  orthonormalization  is  carried  out  over  the  semi- inf inite  interval 
and  is  based  upon  the  Parseval  relation  for  an  inner  product  between  two 


time  functions.  Specifically,  the  Parseval  relation  is 


/  fg*  dt  =  /  F(s)  G*(-s) 

.00  Q 


(34) 
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It  follows  that  the  orthogonality  in  the  time  domain  is  equivalent  to 

orthogonality  in  the  frequency  domain.  If  f  and  g  are  sums  of  complex 

* 

exponentials,  they  will  be  orthogonal  provided  F(s)»  G  (-s)  is  a  rational 

function  which  is  analytic  (i.e.  has  no  poles)  either  in  the  left  half 

or  in  the  right  half  of  the  s-plane.  Consider  the  set  of  exponentials 
s,t 


{e  1  },  i=l ,2, . . . ,n, . . . ,m  t  >  0,  o  <  0 


(35) 


th 


and  construct  an  orthonormal  set  where  the  n  orthonormal  basis 


function  is  given  by 


¥  (s)  = 
n 


s-s 


n-1  s+Sj 

n 


i-i  s"si 


(36) 


th 


Eq.  (36)  can  be  interpreted  in  terms  of  passing  the  n  exponential  time 
function  through  an  all-pass  filter  structured  from  the  previous  n-1 
exponents.  The  all-pass  filter  interpretation  points  out  that  it  is  the 
phase  which  is  responsible  for  orthogonality. 

Thus  far,  orthogonality  has  been  considered  over  a  semi-infinite 
interval.  Provided  a  finite  interval  is  suitably  long,  orthogonality 
can  still  be  approximated  by  this  procedure.  From  (36),  note  that 

s .  t 


y  (t)  =  (s)} 

n  n 

-1 


IV  1 

i=i  1 


(37) 


where  L  is  the  inverse  Laplace  transform. 

In  the  time  domain  orthonormality  requires 


/4'_  (t)H'  (t)  dt 


£=q 

Wq 


(38) 


Substitution  of  (37)  into  (38)  yields 


/V OV  (t)  dt 


£ 

I 

i-1  k-1 


/e(si+sk)  *dt 


i  l -AA: 

i-l  k-l  s1+sk 


£-q 

£j*q 


(39) 
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For  a  finite  interval  of  length  A,  consider 
A  Z  q 


*  <Si+Sk)A  , 

/  V £(t)»  (t)  dt  =  l  l  A 
0  *  <  i=l  k=l  1  k  s  +s. 


(40) 


i  k 


Clearly,  (40)  reduces  to  (39)  provided 
(si+Sk)A 

e  «  1  (41) 

Assume  | a ^ |  <  jo^l  for  i  =  2,3,. ..,n.  Then  (41)  can  be  replaced  by  the 
familiar  inequality 

r  A 

-2y 


(si+sk)A 


2o  A 

<  e  1  =  e 


«  1 


(42) 


Provided  the  interval  length  is  chosen  such  that  (42)  is  satisfied,  the 
orthonormal  basis  generated  by  (36)  is  very  close  to  being  orthonormal 
over  the  finite  interval.  Interestingly  enough,  |y|  £  2,  which  was  the 
condition  for  | pj 2 (A) |  =  |p12(°°)|,  also  satisfies  the  orthogonality 
condition  of  (42) .  Thus 

IyI  1  2  (43) 


is  the  condition  to  be  satisfied  if  a  finite  interval  Is  to  be  considered 
as  though  it  were  an  infinite  interval. 

4.  Conclusion 

_2 

When  |y|  >_  2  or  the  record  length  is  greater  than  —  (where  o  is  the 

°1  1 

real  part  of  the  most  slowly  decaying  exponential)  then  the  finite 
record  length  may  be  considered  as  though  it  were  infinite  and,  in 
addition,  a  set  of  nearly  orthogonal  basis  can  be  generated  over  the 
finite  interval. 
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Fig.^f.  Equal  correlation  magnitude  contours  in  the 

a,  Jg|  plane.  a  >  1  is  assumed  in  the  text. 
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